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ABSTRACT 

Details of various unknown physical processes during the cosmic dawn and the epoch of 
reionization can be extracted from observations of the redshifted 21-cm signal. These 
observations, however, will be affected by the evolution of the signal along the line-of- 
sight which is known as the “light-cone effect”. We model this effect by post-processing 
a dark matter body simulation with an 1-D radiative transfer code. We find that 
the effect is much stronger and dramatic in presence of inhomogeneous heating and 
Lya coupling compared to the case where these processes are not accounted for. One 
finds increase (decrease) in the spherically averaged power spectrum up to a factor of 
3 (0.6) at large scales (fc ~ 0.05 Mpc”^) when the light-cone effect is included, though 
these numbers are highly dependent on the source model. The effect is particularly 
significant near the peak and dip-like features seen in the power spectrum. The peaks 
and dips are suppressed and thus the power spectrum can be smoothed out to a large 
extent if the width of the frequency band used in the experiment is large. We argue 
that it is important to account for the light-cone effect for any 21-cm signal prediction 
during cosmic dawn. 

Key words: intergalactic medium - radiative transfer - cosmology: theory - cosmol¬ 
ogy: dark ages, reionization, first stars - galaxies: formation - X-rays: galaxies 


1 INTRODUCTION 


Understanding the physical processes during the ‘cosmic 
dawn’, when the very first sources formed in the universe, 
and the epoch of reionization (EoR), when the primordial 
neutral hydrogen (H l) got ionized, are some of the major 
goals of present day observational astronomy. The existing 
probes li ke observations of co s mic microwave background 
(CM B l dKomatsu et all I 2 OIII : IPlanck Collaboration et al.l 
[ 201 ^ 201E ) and high redshift quasar ab s orption spectra 
llGunn fc Petersc^ll96d : iBecker et al.ll200il : IfG et ani2003l . 
I 2 OO 6 I: iGoto et al.ll20lif ) suggest that the process of reion- 
izat ion occurred mainly during the redshift range 6 < z < 
15 llMa lhotra fc Rhoad j I 2 OO 6 I: CTioudhur yfc Ferrari |2006| : 


iMitra et al. 2011 . 20121 : Robertson et al. 2015l l. However, 
these observations are still unable to constrain the details 
of the reionization very accurately, e.g., the exact duration 
of reionization, nature of the reionization process, nature of 
sources present during those epochs, ionization and thermal 
states of the intergalactic medium (IGM) etc. 
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Observations of redshifted H l 21-cm signal from these 
epochs promise to provide us with the details of many such 
unknowns (see Fujlanetto_et_ai] I 2 OO 6 I : [Morales fc Wvith3 
I 2 OIOI : [Pritchard fc Lo^ 2012[ . for details). This signal not 
only carries information about the ionization state of the 
hydrogen in the IGM, but also the thermal state of the IGM. 
The first generation low frequ ency radio telescopes, like Low 
Frequency Array (LOFAR|j |[van Haarlem et aI.|[20ljL the 
Precisio n Array for Probing the Epoch of Reionization (PA- 
PER)0 ([Parsons et al.|[2 01^1. the Murchison Widefield Ar¬ 
ray (MWA'FI ( Bowman et al.][2013l : pTingav et ^ 2013lL the 
Giant Metrewave Radio Telescope ('GMRTi ['^[ ([Ghosh et al.[ 
[ 2 OI 2 I: [Paciga et al.][2013ll , have already started observations 
with the aim to detect the 21-cm signal. These instruments 
are focussing mainly on the EoR as they do not contain the 
lower frequency bands required to probe the cosmic dawn. 
The next generation telescopes like the Square Kilometre 
Array (SKAR which are equipped with lower frequency de- 


^ http://www.lofar.org/ 

^ http://eor.berkeley.edu/ 

® http://www.mwatelescope.org/ 
http://www.gmrt.tifr.res.in 
® http://www.skatelescope.org/ 
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tectors and extremely sensitive, should be able to probe the 
cosmic dawn and reveal the unknown properties o f the very 
first sources in the Universe jMellema et al.ir2013ll . 

The 21-cm signal from H l is sensitive to various 
quantities such as the number density and clustering of 
sources, their ionizing, heating and coupling efficiencies, 
the escape fraction of photons at various frequency bands, 
various feedback effects etc. These quantities are eventually 
parametrized in terms of the collapse fraction of dark 
matter haloes, fraction of baryons converted into stars, 
the stellar initial mass function, Ly a, UV and X-ray 
luminosities, X-ray spec t ral in d ex etc jMorales fc Wy ithd 
I 2 OIOI : IPritchard fc Loe9 l2012l : iMesinger et al.l 120141 '). In 
addition, line of sight effects such as the peculiar velocities 
and evolution of the 21-cm signal could also affect the 
signal significantly. It is thus very important to model 
the expected H l 21-cm signal properly in order to in¬ 
terpret the observations. A large amount of theoretical 
modelling and simulations is needed to explore all pos¬ 
sible scenarios and such analyses are expected to play 
a major part in designing 21-cm experiments. Different 
approaches such as analytical (e.g., iFurlanetto etjalj JgOOj: 


20071: 

Mesineer & Furlanettd 20071: Santos et al. 20081: 

Thomas et al.l 20091: IChoudhurv et al.ll2009l:lLa Plante et al. 

2014 


Ghara et al. 2015|L and numerical (llliev et al. 

2 ooe 


iMellema et al.l I 2 OO 6 I: McOuinn et al.l 20071: 


Shin, Trac fc Cenl l2008l : iBaek et al.l 2009l f have been 


attempted in modelling the signal. 

In this paper we investigate one of the important line 
of sight effects i.e, the ‘light-cone effect’ on the 21-cm signal 
from the cosmic dawn and the EoR. While simulating the 
signal, it is generally assumed that every part of a simula¬ 
tion box has the same redshift. We call each such simulation 
snapshot as ‘coeval box’. In reality, regions which are nearer 
to the observer along a line of sight will have lower redshifts 
than the regions which are far away. As a result the ob¬ 
served signal will have effects of redshift evolution imprinted 
on it. Studies through analytical modelling and simulations 
have been done to understand this so called ‘light-cone’ ef¬ 


fect on the two point corr elation function (iBarkana fc Loeb 


200^ Zawadn^^^^ 2014|) an d power spectrum ( Datta^^^^ 
2012h : lLa Plante et al.ll2014l : [Datta et ^120141 '). Datta et al 


||2014i find that, depending on the observational bandwidth, 
the effect could enhance the power spectrum by a factor of 
up to ^ 5 at the initial stages of the EoR and suppress by 
a significant amount at the last stages of the EoR. Interest- 
ingly, no sig:nifican t light-cone anisotropy has been found in 
iDatta et al.l (l2014h . It has also been noticed that the light- 
cone effect on the power spectrum will be significant when it 
evolves non-lin early with redshift as the linear evolution is 
smoothed out llDatta et al.ll2012bl 'l. However, these studies 
assume that the entire IGM is always heated significantly 
higher than the CMB brightness temperature and that the 
Lya coupling is very strong. These assumptions may hold 
during the later stages of reionization but may not be true 
at the cosmic dawn and the initial phase of reionization. 
It is believed that the inhomogeneities in heating and Lya 
coupling will significantly influence the signal during cosmic 
dawn and early phase of the EoR (when the mass averaged 
ionization fraction a;Hii < 0.2). These make the evolution of 
the H I power spectrum with redshift more dramatic and 


thus we expec t a ver y strong light-cone effect on the signal. 
IZawada et al.l (l2014l 'l too have studied this using numerical 
simulations with the main focus on the two point correla¬ 
tion functions. In this work, we self consistently calculate the 
Lya coupling, heating of the IGM for various source models 
and reionization histories and study the light-cone effect on 
the H I power spectrum. 

We have organized the paper in the following way. In 
section [2l we have briefly described the procedure to gen¬ 
erate the brightness temperature maps of the 21-cm signal 
self-consistently using an 1-D radiative transfer code and 
the method to incorporate the light-cone effect in the co¬ 
eval box. We have presented our results in section [3] be¬ 
fore summarizing and discussing our results in section 21 
The cosmological parameters used for the simulation are 
Dni = 0.32, Da = 0.68, Db = 0.049, h = 0.67, = 0.96, 

and (78 = 0.83, consistent with the recent results of Planck 
mission dPlanck Gollaboration et al.ll20lj L 


2 SIMULATIONS 
2.1 The 21-cm signal 

The redshifted 21-cm signal is expected to be measured as an 
offset from the CMB radiation. The differential brightness 
temperature, observed at a fr equency Vohs along a l ine of 
sight n, can be written as (e.g, IFurlanetto et al]l2006l ) 


5Tb(!Zobs,h) = 5Tb(x) = 27 a;Hi(^:,x)[l + 5 b(2,x)] 


1^0.023 


/ 0.15 1-f 



10 

1 

Ts{z,yi)_ 


mK, 


( 1 ) 


where x is the distance along the line of sight to redshift z, 
which can be denoted as l-\-z = 1420 MHz/j/obs. The quanti¬ 
ties xm{z, x) and Sb{z, x) denote the neutral hydrogen frac¬ 
tion and the density contrast in baryons respectively at point 
X and redshift 2 :. The quantity Ts is the spin temperature of 
the neutral hydrogen gas in the IGM, which is determined 
by the coupling of the neutral hydrogen gas w ith the CMB 
phot ons by Thomson scatterin g, Lya coup ling (|W(QuthuYsen 
19M) and collisional coupling |Fieldin~95a : IFurlanetto fc Oh 
200£). The quantity T~f{z) = 2.73 x(l -|- 2 ) K is the CMB 
brightness temperature at redshift z. Note that we have 
not considered any effect from the peculiar vel ocities (e.g, 
iBharadwai fc aH l2004l : IBarkana Loe^ l2005a^ of the gas 
in the IGM in the above equation. 

The first generation radio telescopes are not sufficiently 
sensitive to image the differential brightness temperature 
of the 21-cm radiation. Instead these are expected to mea¬ 
sure the signal in terms of statistical quantities like the 
power spectrum of 5T\^ fluctuations. We have presented 
our results in terms of the dimensionless power spectrum 
A^(fc) = P{k)/2'K^ which also represents the power per 

unit logarithmic interval in scale fc. Here P{k) is the spheri¬ 
cally averaged power spectrum of STh fluctuations, which is 
defined as 

(5rb(k)5tb*(k')) = (27r)®fo(k - k')P(fc), (2) 

where 5Tb (k) is the Fourier transform of 5rb(x) defined 
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in equation ©• In general the universe is expected to be 
isotropic at large scales, though several effects can gener¬ 
ate anisotropy in the observed signal. Line of sight pecu¬ 
liar velocity of gas (or the so called redshift-space distortion 
(RSD)) in the IGM can affect the 5Th fluctuations and hence 
the signal along the line of sight and thus can introduce a 
difference in the power s pectrum along and perpendicular 
to th e observed direction (iMao et al.ll2012l : lMaiumdar et aP 
I2OI3I : I Jensen et al.l[ 201 ^ . On the other hand the Alcock- 
Paczynski effect can also signif icantly contribute to the 
anisotropy in the observed signal (lAjcock_fc Paczynskill 19791 : 
lAli, Bharadwai fc Pand^ l2005l : BarkanaT 2006l i. In gen¬ 
eral, the anisotropic 21 -cm power spectrum is denoted as 
A^(fc,/r), where fi = cos 9 with 9 being the angle between 
the line of sight and the Fourier mode k. Since the light-cone 
effect too influences the signal only along the line of sight, 
it can in principle make the signal anisotropic. In this study 
we mainly focus on the light-cone effect on the spherically 
averaged power spectrum and we will discuss the anisotropy 
introduced by this to the 21 -cm power spectrum in section 


The simulation of the STh signal involves three main 
steps: (i) generation of underlying baryonic density and ve¬ 
locity fields, (ii) modelling the radiation sources and (iii) 
computing the propagation of ionization and heating fronts 
using an 1-D radiation transfer code. The method u sed in 
this paper is essentially described in iGhara et al. I ll2015l. 
hereafter Paper I). We summarize the main steps in the 
next few subsections. 


2.2 A—body simulations 

The density and velocity fields, used in our simulatiom 
are generated using the pub licly available code CUBEP^^ij 
dHarnois-Deraps et al.ll2013l ~) which is a massively parallel 
hybrid particle-particle-particle-mesh (P^M) code. The sim¬ 
ulation started at redshift z = 200 with the initialization 
of the particl e positions and ve locities using CAMB trans¬ 
fer functiorQ (iLewis et al. I l 2000 l i and employing Zel’dovich 
approximation. 

The properties of the simulations we have used in this 
work are: (i) the number of particles used is 1728®, (ii) the 
box size is 200 h~^ comoving Mpc (cMpc), (iii) the num¬ 
ber of grid points used are 3456® and (iv) the mass of each 
dark matter particle is 2 x 10® M©. Since the position and 
velocity arrays of the simulation are too large in size, we 
prefer to generate the density and velocity fields on a grid 8 
times coarser than the simulation grid using a top-hat filter. 
These snapshot files are generated between 25 ^ 2 > 6 in 
equal time gap of 10^ years. The baryons are assumed to be 
simply tracing the dark matter, i.e., if dark matter density 
at position x is pdm(x), then the baryonic density will be 
Pb(x) = (fiB/Um) X pdm(x). The velocity files on the coarse 
grid are used to incorporate the peculiar velocity effects in 
the 21 -cm brightness temperature. 

For each snapshot, haloes were identified using a run¬ 
time halo finder algorithm which uses spherical over-density 
method. The minimum halo mass resolved using this method 


is 2 X 10® M 0 . Apart from these haloes, it is expected 
there will also be considerable number of small mass haloes 
^ 10® Mq where the gas can cool via atomic transitions and 
form stars. As the number of such small mass haloes is quite 
high, these haloes may have significant impact on the reion¬ 
ization scenarios, particularly at the early stages. Identify¬ 
ing such small mass haloes with the spherical over-density 
method requires simulation box of a higher resolution, which 
is somewhat beyond our ability right now. Hence we employ 
a sub-grid method to find the haloes down to masses as 
small as 10® Mq . We have followed the extended Press- 
Schech ter model of BondetjalJ (Il99ll i and hybrid prescrip¬ 


tion of iBarkana fc LoebI ( 2004') for implementing the sub¬ 
grid model (for details please see Paper I). 


2.3 Modelling the sources 

The main sources of ionizing photons are the stars residing in 
galaxies, which form within dark matter haloes. We assume 
that haloes having mass > 10® Mq contribute to ionizing 
radiation in regions that are neutral. Within ionized regions, 
the corresponding threshold mass is taken to be 10® Mq 
because of radiative feedback. The stellar mass of a galaxy 
in a hosting dark matter halo of mass Mhaio is 

M. = L Mhaio, (3) 

where /* is the fraction of baryons residing within the stars 
in a galaxy. This fraction depends on the metallicity and 
mass of the galaxy and is not well constrained for early 
galaxies. Hence for simplicity we have assumed a constant /* 
throughout the reionization history, with its value is chosen 
in such a way that the resulting reionization history is con¬ 
sistent with the constraints obtained from the optical depth 
measurement of the CMB observations. 

The spectral energy distribution (SED) of a galaxy can 
be calculated using stellar population synthesis codes pro¬ 
vided one knows the values of the initial mass function 
(IMF) of the stars, initial metallicity etc. We assume a 
Salpeter IMF for the stars with mass range 1 to 100 Mq 
throughout the reionization histo ry. The meta l licity e volu¬ 
tion is taken from the mode ls of IPaval et al.l (l2009al 'l and 
IPavak Ferrara fc Sard ll 2010 ll which are consistent with the 
best-fit mass-metallicity relation. As in Paper I, we gener¬ 
ate the UV and near-infrared (NIR) spectral energy distri¬ 
butions of the galaxies for standard star forma tion scenarios 
using the code PEGASEifl dFioc fc Roccalfl997lf . The fraction 
of UV ionizing photons that escape into the IGM from the 
galaxies (/esc) is highly uncertain for high redshift galaxies. 
For every redshift of interest, we choose the value of /esc so 
as to obtain a given reionization history, which is explained 
later in section [3] 

In this simulation, we have considered ionization of both 
hydrogen and helium in the IGM by central sources. The 
photons from the stars can only ionize H l and He l, however 
ionization of He ll and heating of the IGM is not efficient 
by the stellar sources because of the lack of high energy (> 
50 eV) photons. Following our previous work (Paper I), we 


° http://wiki.cita.utoronto.ca/mediawiki/index.php/CubePM 
^ http://cainb.info/ 


http://www2.iap.fr/pegase/ 
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model the X-rays from galaxies as having a power-law SED 
with a spectral index a: 

I,{E) = A (4) 


where the normalisation constant A can be fixed in terms 
of the ratio fx of X-ray to UV luminosity from the 
source. This kind of a X-ray spectrum is expected from 
a miniquasar-like source, i.e., a galaxy with intermediate 
mass black hole of mass in the range of 10^ — 10^ Mp 


L U iVip 

i]l200ll: 


llElvis et al.lll993:lLaor et al.lll997l:IVanden Berk et al.ll 

IVignali, Brandt, fc Schneider! l2003[ l. There could be other 

sources of high energy phot on s like the high mass X-ray bi¬ 
naries jFragos et al.l I^13al f9: iFialkov et al.l [20lJ}_ and the 
hot in terstellar medium within early galaxies ( Pacucci et al.l 
l2014l ~l. However, we concentrate on a power-law X-ray spec¬ 
trum in this work. 

Note that the parameter fx can be related to the cen¬ 
tral black hole to galaxy mass ratio. Since main purpose of 
this paper is to study the light-cone effect for models which 
incorporate heating and Lya coupling self-consistently, we 
have kept the values fx and a to be fixed. We choose fx to 
be 0.05 and a to be 1.5 for our study, which are consistent 
with the observations of local quasars. The effect of choosing 
other values of fx and a was studied in detail in Paper I. 


2.4 STh maps using 1-D radiative transfer 

As soon as the first sources of light appear in the dark matter 
haloes, they start to influence the ionization and the ther¬ 
mal state of the IGM in the universe. Since the UV photons 
have relatively smaller mean free path, they will be absorbed 
by the neutral hydrogen in the immediate neighbourhood, 
resulting in highly ionized bubbles around the sources. How¬ 
ever, the X-rays can propagate long distances in the IGM 
and partially ionize and heat up the medium. Although it is 
necessary to carry out 3-D radiative transfer simulations to 
generate the ionization and heating maps, it requires huge 
computational power and long run time. We rather prefer an 
alternative method based on 1-D radiative transfer which is 
faster and hence more efficient for expiring the unknown 
param eter space (iThomas et al.l 120091 : iThomas &: Zaroubil 
I 2 OIII I. The method is briefly described in the following, we 
refer the readers to Paper I for the details. 

(i) The ionization and heating profile around a galaxy 
depends on the luminosity and the surrounding neutral hy¬ 
drogen distribution. We generate the ionization and heating 
profiles around sources for different galaxy masses, redshifts 
and background gas densities. Though the baryonic distri¬ 
bution around a galaxy changes with distance, but for sim¬ 
plicity we assume the gas is distributed uniformly around 
the source. This is one of the limitations of our method. 

(ii) Given the list of haloes at a certain redshift, we cal¬ 
culate the corresponding galaxy masses within the haloes 
as stated in section We assume that no star formation 
occurs within low mass haloes (< 10® Mq). In order to in¬ 
corporate radiative feedback, we suppress the galaxy forma¬ 
tion within newly formed haloes having mass M < 10® Mq 
if they form in already ionized regions, i.e., regions of the 
IGM with xmi larger than 0.5. 0 

® The radiative feedback is effective when the IGM is photo- 


(iii) Depending on the mass, background gas distribu¬ 
tion and redshift of the galaxies, we estimate the ionizing 
photons from the galaxies. The ionized bubbles around indi¬ 
vidual sources are computed by solving the radiative transfer 
equations in 1-D. In case of overlap between the bubbles, we 
estimate the unused ionizing photons and distribute them 
equally among the overlapping ionized bubbles. 

(iv) The X-rays can penetrate into and partially ionize 
the neutral IGM. Beyond the highly ionized regions around 
the sources, the ionization fraction is estimated using the 
pre-generated profiles. Once ionization maps are generated, 
the kinetic temperature maps are generated using a corre¬ 
lation of kinetic temperature and ionization fraction in the 
partially ionized regions (see Paper I for details). 

(v) As far as the radiative transfer of the Lya pho¬ 
tons are concerned, we simply assume the escape fraction 
of the Lya photons to be 1 and the number density of Lya 
photons from the galaxies decrease as 1/df, where d be the 
distance from the source. The Lya coupling coefficient is cal¬ 
culated from the Lya flux, which completes every quantity 
required for computing the Tg maps. The brightness tem¬ 
perature maps are generated in the simulation box using 
equation O- 

(vi) The effect of the peculiar velocity of the gas in the 
IGM is incorporated into the STh maps using the cell move¬ 
ment me thod (or M esh-to-Mesh (MM)-RRM scheme) de¬ 
scribed in iMao et al.l (l2012l L which is time efficient and com¬ 
putationally easy to implement. The details can be found in 
Paper 1. 


2.5 Implementing the light-cone effect 

The light-cone effect essentially accounts for the evolution 
of the signal with redshift wit hin the simu l ation b ox or the 
so called “coeval cube” fGGL lPatta et aH (l2012lJ l have de¬ 
scribed the procedure to incorporate the effect and produce 
the light-cone cube. In our study we have incorporated the 
effect in individual coeval cubes instead of ge nerating the 
full red shift range light-cone cube as done in IPatta et al.l 
l|2012blL We assume that the central “slice” of the coeval 
cube represents the cosmological redshift, where a “slice” is 
a 2-D map having width ~ 0.7 cMpc (corresponding to one 
pixel of the simulation box). Along a given line of sight, dif¬ 
ferent distances represent different redshifts, and thus STh 
will be different from that in the coeval cube. The steps to 
include the light-cone effect are as follows: 

(i) We generate the so called “coeval cubes” of 21-cm 
STh maps at different redshifts following the method stated 
in section Let us assume that we have N numbers of 
such cubes which correspond to the cosmological redshifts 
Zi, where i = 1,N, with zi{zn) being the lowest (highest) 


heated to a temperature > lO'^ K, which raises the cosmological 
Jeans mass and thus suppresses galaxy formation in low mass 
haloes with mass M < 10^ Mq. In general, regions with temper¬ 
ature > 10"^ K are highly ionized (with xhii > 0.5), hence our 
implementation of feedback should be reasonable. In more realis¬ 
tic scenarios, the feedback is not expected to be so drastic, e.g., 
M < 10®Mq haloes within ionized regions can retain a fraction 
of their gas and possibly form stars at a lower rate. Our method 
does not account for such effects. 
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redshift in our simulation (in our case, we have 77 such co¬ 
eval cubes between redshift 6 and 20). 

(ii) Each of these coeval cubes have a length L and 
contain rd grid points, which in our case have the values 
200 h~^ cMpc and 432^ respectively. We assume the central 
2-D slice of the cube, in the perpendicular plane of the line 
of sight (i-axis in this case), correspond to the cosmological 
redshift of the coeval cube. The comoving distance D{zi, zL) 
between two redshifts zi and Zu is given by. 







where c is the speed of light and H{z) is the Hubble con¬ 
stant at redshift z. Depending on the distance from the 
central slice along the line of sight, we calculate the cor¬ 
responding redshift for each slice. For example, the redshift 
Zp corresponding to the pth slice is calculated by demand¬ 
ing that D{zc,Zp) in equation ((h]) is equal to the distance 
(p — n/2) X Ljn of the pth slice from the central slice. Note 
that Zc is the redshift corresponding to the central slice. In 
the light-cone cube, the signal at each slice would be at the 
redshift corresponding to that slice and not the cosmological 
redshift of the coeval cube. 

(hi) If the pth slice corresponds to a redshift Zp, which 
satisfies the condition Zi < Zp < Zi+i, then the 5Tb maps 
for that slice is computed by linearly interpolating the maps 
of the slices at redshifts Zi and Zi+i. Applying the same 
procedure to all the slices, the “light-cone cube” is generated 
which corresponds to the redshift of the central slice. 

Note that the redshift-space distortion needs to be ap¬ 
plied to the brightness temperature maps before implement¬ 
ing the light-cone effect, otherwise the redshift evolution of 
the gas velocities will not be incorporated properly. 


3 RESULTS 

In order to compute the 5Tb maps, one Erst has to choose 
appropriate reionization and heating histories, which in turn 
depend on the nature of sources. In this work, we consider 
three models of radiation sources S'!, S2 and S3 as listed in 
Table [T] to investigate the eSect of difierent reionization his¬ 
tories on the 5Tb maps. For all the three source models, we 
assume each galaxy to contain a mini-quasar like source at 
the centre which radiates X-rays with a power law distribu¬ 
tion having index a =1.5, and the X-ray to UV luminosity 
fraction of the source fx is 0.05. We have fixed the pa¬ 
rameter /» to 0.1 throughout the reionization history, while 
the escape fraction /esc is taken to be different in different 
models. In model 51, we take /esc = 0.1 which causes the 
reionization to start at 2 : ~ 20 and complete at 2 : ~ 6.5. 
The electron scattering optical depth r turns out to be 0.07, 
which is consistent with the recen t constraints from Planck 
dPlanck Collaboration et al.ll2015l l. The mass averaged neu¬ 
tral hydrogen fraction, shown in Figure [T] attains a value 
0.99 and 0.5 around redshift 14 and 8.5 respectively for 51. 

The value of /esc = 0.2 is higher in model 52 lead¬ 
ing to early reionization completing at 2 : ~ 8. The result¬ 
ing value of T = 0.076 is higher than that in model 51 
but is still consistent with the recent Planck constraints 



Figure 1. The solid curves represent the evolution of the mass 
averaged neutral fraction of hydrogen (xjji) in the IGM. The 
dashed curves represent the evolution of (1 — Cheated), where 
Xbeated is the volume averaged heated fraction of the IGM. The 
thin and thick curves represent models 51 and 52 respectively. 
We have defined ‘heated region’ as a region with Tk larger than 
Try. The X-ray spectrum from the central mini-quasar follows a 
power law with index o =1.5 and the X-ray to UV luminosity of 
the source is 0.05 for this model. 


llPlanck Collaboration et al.l I 2 OI 5 I I . In presence of a mini¬ 
quasar like X-ray source at the centre of each galaxy, the 
evolution of the heated fraction (defined as the volume frac¬ 
tion of regions with Tk > Try) is much faster than the evolu¬ 
tion of ionization fraction. As shown in Figure[T] the universe 
becomes completely heated below redshift ~ 12 for both the 
source models 51 and 52. 

Note that the minimum halo mass for 51 and 52 source 
models is ^ 2 x 10® Mq, i.e., we have not included any 
sub-grid sources in these two models. The effect of small 
mass haloes is studied in model 53 where haloes as small 
as 10® Mq are included using the sub-grid prescription dis¬ 
cussed in The value of /esc is varied at every redshift 
so as to generate a reionization history (i.e., mass averaged 
neutral fraction) identical to the model 52. In addition, we 
have included the effects of radiative feedback in model 53. 

As in paper I, we will consider three different models of 
heating and Lya coupling and study the light-cone effect. 
The models are listed in Table [21 Model A assumes Ts to be 
m uch larger than Try, t his model is identical to models used 


Datta et al. l (|2012bD ;[Lr Plante et al.l ll2014l '): [Datta et all 
(l2014h . Model B assumes the IGM kinetic temperature to 
be coupled with the spin temperature through Lya radi¬ 
ation (i.e, the Lya coupling coefficient Xa{x) 1), but 
the heating is calculated self-consistently. Finally in model 
C, we have calculated the heating and Lya coupling self- 
consistently and thus it is the most realistic model to be 
considered while generating 5Tb maps. The ionization his¬ 
tory in Figure [T] is same for all the three models A, B and 
C, while the heating history is same for the models B and 
C. The effect of peculiar velocities of the gas is also taken 
into account in all the three models. 


© ? RAS, MNRAS 000. nUT^ 





















6 Ghara, Datta & Choudhury 


150 

j- 100 

>. 

50 



0 50 100 150 200 

X (cMpc/h) 


0 50 100 150 200 

X (cMpc/h) 


200 


Figure 2. Randomly chosen two-dimensional slice of the 6T\^ distribution with RSD included from the coeval cube (left) and light-cone 
cube (right) at redshift 9.5 (with mass averaged ionization fraction 0.5) for the model A with the source model S2. The comoving size of 
the box is 200 cMpc h~^. The central redshift of the light-cone (LC) cube is 9.5, whereas the redshift span from left to right along the 
x-axis (LOS direction) is 8.86 to 10.13. Note that this model makes the assumption Tg >> Ty and is driven by haloes identified using 
spherical overdensity halo-finder in the simulation box. 
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Figure 3. Similar to Figure[2] but for the model C and the coeval cube corresponding to redshift 13 when the universe was 50 % heated 
by volume for the source model S2. The left and right hand edge of the ‘LC’ box correspond to redshifts 12 and 14 respectively. 


3.1 Visualizing the light-cone effect 

The light-cone effect can be easily visualized in the maps of 
5Th distribution for different models. Figure [2] shows the 5Th 
map computed on a randomly chosen slice from the simula¬ 
tion box for the model A with the source model S2. The left 
hand panel shows the distribution of 5Tb from the coeval 
cube (with redshift space distortion included, but without 


the light-cone effect) at redshift 9.5 with mass averaged ion¬ 
ization fraction 0.5. The right hand panel shows the same 
slice with the light-cone effect incorporated. Along the line of 
sight (horizontal axis), the central redshift of the light-cone 
(LC) cube is 9.5, whereas the light-cone cube spans from 
redshift 8.86 (left edge of the box) to 10.13 (right edge of 
the box). The ionization bubbles (black regions) are larger 
at the left-hand side of the middle of light-cone cube as 
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Figure 4. Model A: The light-cone effect on the power spectrum of 21-cm brightness temperature fluctuation. The upper most panels 
show the spherically averaged power spectrum of 5T\^ with (thin curves) and without (thick curves) light-cone effect as a function of 
redshift for four different scales fc/Mpc~^ = 0.05 (solid), 0.1 (dotted), 0.5 (dashed) and 1.0 (dot-dashed ). The middle row panels show 
the ratio of the power spectra with and without light-cone effect as a function of redshift. The bottom most panels show the ratio of 
the power spectra with and without light-cone effect as a function of scales for different ionization fractions (0.1,0.3, 0.5, 0.7 and 0.8). 
This model makes the assumption that Tg >> T.y. Three columns represent three different source models 51,52 and 53 respectively. 
Redshift space distortion is included in all the panels. 


compared to the coeval cube. This is because of the fact 
that the left-hand side regions correspond to lower redshifts 
than the middle region of the box and thus correspond to 
later stages of reionization where the bubbles are of system¬ 
atically larger. It is opposite in the case of right-hand side 
regions. 

Figure [3] shows similar slices at redshift 13 for the model 


C. In this case not only the ionized regions appear larger in 
the front side of the light-cone cube compared to the co¬ 
eval cube, but the heated regions (regions with Tk > T^) 
too show similar trend. One should notice that the left 
half shows the signal predominantly in emission whereas the 
right half shows the same mostly In absorption, although the 
corresponding coeval slice shows that emission and absorp- 
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Source 

model 

/esc 

r 


M^in 

halo 

(Mo) 

sub-grid 

halo 

radiative 

feedback 

51 

0.1 

0.07 

6.5 

2 X 10® 

no 

no 

52 

0.2 

0.076 

7.8 

2 X 10® 

no 

no 

53 

varying 

0.076 

7.8 

10® 

yes 

yes 


Table 1. Properties of the source models considered in this paper. 
The UV escape fraction /esc for the model 53 is varied to get 
similar ionization history like 52. 


Model 

Lya coupling 

Heating 

A 

Coupled 

Heated 

B 

Coupled 

Self-consistent 

C 

Self-consistent 

Self-consistent 


Table 2. Properties of models considered to investigate the effects 
of inhomogeneous heating and Lya coupling on the signal. The 
terms ‘coupled’ and ‘heated’ represent the scenarios where Xa 3> 
1 and Tk ^ Tj respectively. 

tion regions are homogeneously distributed over the entire 
slice. This significant imbalance of the emission regions be¬ 
tween the left and right-hand regions of the light-cone cube 
is expected to introduce effects on the observable signal at 
initial stages of reionization. 

3.2 Comparison with previous studies 

It is possible to validate our formalism by comparing the 
model A where we have assumed Ts » with the exist¬ 
ing works in the context of the light-cone effect on the 21-cm 
signal. Figure [4] shows the light-cone effect on the spheri¬ 
cally averaged power spectrum of the 21-cm 5Ti, fluctuations 
(with redshift space distortion included) for the model A0 
The upper panels show the evolution of the dimensionless 
power spectrum for the three source models as a function 
of redshift for four different scales k = 0.05, 0.1, 0.5 and 1 
Mpc“^ which are represented by the solid, dotted, dashed 
and dot-dashed curves, respectively. The signal is always in 
emission for this model and thus the power spectra have rel¬ 
atively smaller amplitudes 10 mK^ ) at large scales. The 
large scale power spectrum for each of the three source mod¬ 
els shows one trough and one peak when plotted as a func¬ 
tion of redshift around the period when universe was ~ 10% 
and ^ 65% ionized by mass. Initially the highly dense re¬ 
gions get ionized and thus H l fluctuations at large scales 
are suppressed. This results in the tr ough like feature in the 
evolution of the p ower spectrum fe.g. lMaiumdar et al.ll20i^ : 
iDatta et al.ll2014l '). The power spectrum then increases with 
the growth of bubble size and finally decreases with the de¬ 
crease in ionization fraction, which generates the prominent 
peak around ionization fraction 0.65 for this model. 

Note that we exclude modes with ky = 0 while calculating 
the power spectrum, as these modes are not accessible to the 
interferometers. 


The thick (thin) curves in the upper panels of the fig¬ 
ure represent the power spectra without (with) the light- 
cone effect included. It is clear that the light-cone effect is 
most substantial around the trough and the peak regions in 
the evolution of the power spectra. Because of this effect, 
the power spectra can be enhanced by a factor as large as 
~ 50% (when xmi ~ 0.15) and suppressed by ~ 20% (when 
xhii ~ 0.8) for the model 52 at large scales k = 0.05 Mpc“^, 
as shown in the middle panels of Figure[4l The effect is min¬ 
imum around the period when universe is ~ 50% ionized. 
This is because of the fact that around 50% ionization, any 
linear evolution of the power spectru m as a function of z i s 
smoothed out, as was pointed out bv iDatta et al.l (l2012blf . 
On the other hand, when a;Hii ~ 0.1 and 0.8 the evolution 
of the power s pectrum is highly non-line ar which makes the 
effects strong jPatta et al.ll2012bl . 12013 1. We also see that 
the effects are less for a smoother ionization model like 51 
compared to 52. The reason is that in case of 51, the evolu¬ 
tion of the ionized bubbles are less rapid compared to model 
52. The evolution of power spectrum is more non-linear at 
the trough and the peak regions in presence of small mass 
haloes as shown in the upper panel for the model 53. Thus 
the light-cone effect in model 53 is larger compared to the 
other two models. The bottom most panels of the figure 
show the scale dependence of the ratio of the power spec¬ 
trum with and without light-cone effect at different stages 
of the reionization history. The evolution of the power spec¬ 
trum at small scales is much more linear as shown in the 
upper panels of the figure and thus the light-cone effect is 
smaller at small scales compared to larger scales. Our re¬ 
sults for the model A are cons istent with a similar model 
presented in lPatta et al.l (I20l4 l. 

3.3 Effect of inhomogeneous heating and Lya 
coupling on the light-cone effect 

We now discuss the light-cone effect in presence of non- 
uniform heating and Lya coupling. Figure [S] shows the same 
quantities plotted in Figure U but for the model B (where 
the fluctuations in heating are calculated self-consistently, 
but the Lya coupling process is assumed to be very effi¬ 
cient all throughout). The top panels show the evolution of 
power spectra for different length scales. One can compare 
the plots with those in Figure 3] and immediately conclude 
that the evolution is much more rapid in the model B. For 
example, the amplitude of the power spectrum is ~ 50 mK^ 
at z 20 which increases up to 500 mK^ at z: ~ 15 and 
then decreases rapidly to ~ 0.5 mK^ at 2 ~ 11 for k = 0.05 
Mpc“^ (in the 52 model). The corresponding evolution is 
much less rapid for the model A where the amplitude is al¬ 
most constant at ~ 1 mK^ from z = 20 till 2 ~ 14. The 
light-cone effect is thus expected to be more significant for 
the model B. 

The evolution of the power spectrum is faster in model 
B, due to the fact that the heated bubbles are larger and 
grow much faster than the ionized bubbles. This causes the 
evolution of the heated fraction to be very rapid compared 
to evolution of the ionization fraction (see Figure [T]). As a 
large fraction of the gas in the IGM has temperature Tk less 
than T.y during cosmic dawn and initial phase of reioniza¬ 
tion (as gas temperature falls as (1 -I- z)^ due to adiabatic 
expansion of the universe after the decoupling around red- 
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Figure 5. Same as Figure [J] but for the model B where the IGM is assumed to be Lyo coupled while the heating is calculated 
self-consistently. 


shift Zdec ^ 150 and falls as (1 + z)) and the model 
assumes Ts = Tk, the model shows very large absorption 
signal from the cosmic dawn until the IGM is sufficiently 
heated to follow the model A from 2 ~ 10 for S2. The inho¬ 
mogeneous X-ray heating from the mini-quasars results in 
increasing the 5Th fluctuations, which produces a distinct 
peak in the power spectrum at large scales around 2 ~ 15 
when plotted against 2 . The second prominent peak occurs 
at 2 ~ 9 because of ionization fluctuations analogous to 
the model A. The amplitude of the power spectrum at the 
heating peak 10^ mK^) is significantly higher compared 


to the ionization peak (~ 10 mK^) and hence this redshift 
range might be of interest to interferometers like the SKA. 
The reason for this is that at the ionization peak around 
2^9, the power spectrum is dominated by ionization fluc¬ 
tuations and the signal is in emission because the heating is 
already substantial. On the other hand, at the heating peak 
around 2 ~ 15, the 5Tb held consists of both emission and 
absorption regions, and the power spectrum is dominated 
by Ts fluctuations (Paper I). 

As can be seen from Fig[5l the light-cone effect is much 
more prominent in the model B compared to the model A. 
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Figure 6. Same as Figure (5] but for the model C where both Lyo coupling and heating are calculated self-consistently. 


When the light-cone effect is included in the analysis, we find 
maximum suppression of power spectrum around the two 
peak regions and increase around the trough region between 
the two peaks. Inhomogeneous heating makes the evolution 
of the power spectrum much more non-linear and results in 
a stronger light-cone effect on the power spectrum than that 
in the model A. The figure also shows that the light-cone ef¬ 
fect can increase the power spectrum by a factor as large as 
1.23 and 2.84 times around redshifts 16 and 11 respectively 
at large scales {k ^ 0.05 Mpc“^) for the model S2. The sup¬ 
pression of the power spectrum around the peaks around 
redshifts 8 and 15 can be ~ 20 — 25% at large scales. The 


redshift evolution of the power spectrum at small scales is 
almost linear during the cosmic dawn and initial phase of 
reionization as shown in the upper most panels of the figure 
and thus the light-cone effect is minimum at small scales 
during these epochs. Around redshift 12 when the IGM be¬ 
comes highly heated, the evolution in the small scale (5Tb 
fluctuations becomes extremely fast which results in the in¬ 
crement of the power spectrum by a factor of 1.5 due to the 
light-cone effect at small scales (fc ~ lMpc~^). The light- 
cone effect is less in the model 51, as the evolution of the 
ionized as well as the heated regions is less rapid compared 
to that in the model 52. In presence of small mass haloes 
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Figure 7. The redshift evolution of the difference between the power spectrum with and without light-cone effect at scales A;/Mpc“^ = 
0.05 (solid), 0.1 (dotted), 0.5 (dashed) and 1.0 (dot-dashed) for the model C. Left to right panels represent source models 51,52 and 53 
respectively. 


in the model S3, the heating is delayed compared to the 
model S2 as we have varied /esc and kept fx fixed. In this 
case, the fluctuations in STh at large scales is smaller com¬ 
pared to that in the model S2 and thus the amplitude of 
the heating peak is smaller. As the heated regions around 
the sources are very small in S3 and eventually increase, 
the power spectrum shows a dip in its evolution with the 
redshift as shown in the top-right panel of the figure. The 
light-cone effect is able to increase the power spectrum ~ 2 
times around the dip (at 2 : 14). 

The amplitude and overall nature of the heating peak 
around redshift 15 depend on the nature of X-ray sources. As 
for example, the peak may not present if heating is driven 
by X-ray sources like high mass X-ray binaries (HMXBs) 
as they do not have large number of soft X-ray photons 
and thus the heat ing will be almost homogeneous (e.g, 
iFialkov et al.|[20l3 ). In such cases, STh fluctuations will be 
dominated by the neutral hydrogen fluctuations if we con¬ 
sider the Lya coupling to be very high from the beginning. 
Thus the light-cone effect in such scenarios will be similar 
to the model A. Whereas, for X-rays, from sources like su¬ 
pernova, hot-interstellar gas, follow similar profile like mini¬ 
quasar and thus expected to have similar light-cone effect. 

The difference between models B and C is that the Lya 
coupling is calculated self-consistently in model C whereas 
in model B it is assumed that the IGM is Lya coupled. 
As the condition Xa{'x.) S> 1 is not valid for major part of 
the IGM during the cosmic dawn, the amplitude of 5Tb (x) 
in the model C will be much less than the predicted STh 
in the model B. As a result the amplitude of the power 
spectrum is much lower than in the model B initially. The 
inhomogeneous Lya coupling during the cosmic dawn in the 
model C increases the fluctuations in STh and results in a 
distinct peak in the evolution plot of the large scale power 
spectrum as a function of redshift as shown in the upper 
panels of Figure [S] We note that a very small amount of 
Lya photons is sufficient to couple Ts with Tk of the IGM, 


thus the model G follows the model B very soon the first 
sources formed (around redshift 13). 

In presence of inhomogeneous Lya coupling the evolu¬ 
tion of the power spectrum is more dramatic at both large 
and small scales at cosmic dawn. The non-linear rise of the 
power spectrum, from 1 mK^ to ~ 100 mK^ within red¬ 
shift interval 20 to 16, results in a stronger light-cone ef¬ 
fect at the very beginning of the reionization epoch. This 
increases the power spectrum by a factor of ~ 1.5 around 
redshift 18.5 and suppresses by a factor of 0.6 around red¬ 
shift 16.5 (which corresponds the Lya peak for the source 
model 52). The light-cone effect further enhances the power 
spectrum by a factor of ^ 2 around redshift 15.5 (trough 
region between the Lya and heating peaks), followed by a 
suppression by a factor of 0.7 around redshift 14. For red¬ 
shift 2 < 14, the light-cone effect is similar to the model B, 
only difference being that in presence of inhomogeneous Lya 
coupling the effect is little weaker than that predicted by the 
model B. The power spectrum evolves rapidly at small scales 
at the cosmic dawn, which results in a large light-cone ef¬ 
fect at small scales (^3 — 4 times enhanced). As expected, 
the light-cone effect is found to be smaller in model 51 com¬ 
pared to model 52. The enhancement in the power spectrum 
due to the light-cone effect is ~ 3 around redshift 14 which 
corresponds to the trough region between the Lya coupling 
and heating regions for the model 53. 

Importantly, we also notice that the peaks and dips 
found in the evolution of the power spectra for the mod¬ 
els B and G are smoothed out to some extent due to the 
light-cone effect (upper panels of Figs [5] and [5]). The effect 
is more prominent at large scales fc < 0.1 Mpc“^. The light- 
cone effect lowers the peak height and raises the dip by some 
amount. This is particularly important for 53 in model C, 
which is probably the most realistic reioniz ation model we 
consi der here. It has been suggested (e.g., iMesinger et al.l 
l2014h that the peaks can be used to extract source proper¬ 
ties, put constrain on the X-ray and Lya background during 
the cosmic dawn. We argue that the light-cone effect should 
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Figure 8. The dependence of the light-cone effect on the simulation box size. The left hand panel shows the evolution of the power 
spectrum at a scale fc = 0.1 Mpc“^ for the two different boxes Boxi (200 h~^ cMpc, blue-dotted curves) and Boxs (100 h~^ cMpc, red-solid 
curves) respectively with (thin) and without (thick) the light-cone effect. The right hand panel shows the redshift evolution of the ratio 
of the power spectra with and without the light-cone effect. The plots are for 53 model and redshift space distortion effect is included 
while calculating the power spectrum. 


be considered while extracting those parameters from the 
peak and dip measurements. 

Finally, in Figure [7] we show the redshift evolution of 
the difference in the power spectra with and without the 
light-cone effect for the model C for the three source models 
at different scales. As expected, the difference is maximum 
at the peak and the trough regions found in the evolution 
of the power spectra. At large scales (fc ~ 0.05 Mpc“^), the 
difference can be as large as ~ —100 to 100 mK^ . The 
difference is ~ —250 to 100 mK^ at intermediate scales like 
fc ^ 1 Mpc“^. In principle, such a strong effect should easily 
be detected by future experiments like the SKA. 


3.4 Effect of simulation box size 

It has been found earlier by iDatta et al.l (l2014l ) that the 
light-cone effect is larger when larger simulation box is used. 
In order to quantify the effect of box size on our results, we 
estimated the signal using a smaller box of length 100 
cMpc (in addition to our default box of 200 h~^ cMpc). 
The results are shown in Fig. [H] In the left panel, we plot 
the evolution of at a scale fc = 0.1 Mpc“^ for the two 
boxes with and without the light-cone effect. As expected, 
the power spectra without the light-cone effect (thick lines) 
for the two boxes agree with each other. However, we can 
see that the light-cone effect is conside rably less prominen t 
for the smaller box in agreement with iDatta et al.l (|2014|) . 
In particular, the smoothing of the three peak nature in 
the evolution plot of the power spectrum is more prominent 
for the larger box. The effect can be much m ore stronger 
in simulation box with siz e ~ 600 cMpc (e.g., IDatta et al.l 
I 2 OI 4 I : iMesinger et al.ll2014l ) , which may completely smooth 
out the three peak nature of the evolution plot of the power 
spectrum. This will constrain us to choose smaller frequency 
bands during 21-cm experiments to avoid strong light-cone 
effects and restore the peakiness feature of the power spec¬ 
trum which is very useful for parameter estimation etc. The 
same conclusions can be drawn from the right hand panel 
too where we have shown the ratio of the power spectra 
with and without the light-cone effect for fc = 0.1 Mpc“^. 


Clearly, the ratio deviates from unity quite prominently for 
the larger box compared to the smaller one. 


3.5 Anisotropy from the light-cone effect 


Since the light-cone effect modifies the 21-cm signal along 
the line of sight direction similar to the RSD, it is expected 
that it may cause anisotropies in the signal. This was inves¬ 
tigated in detail bv IPatta et al.l (l2014l ) for a model similar to 
our model A, and they concluded that the light-cone effect 
does not induce any significant anisotropies in the signal. 
We confirm their findings for the model A. In addition, we 
find that the light-cone effect does not cause any promi¬ 
nent anisotropies for the models B and C too for relevant 
scales of interest. In fact, for large scales fc < 0.1 Mpc“^, our 
simulation box does not contain sufficient number of modes 
leading to large sample variance. Hence it is difficult to draw 
any significant conclusion on anisotropies at large scales. 

In order to explain the effect in a somewhat simpler 
man ner, we have calculat ed the ‘anisotropy ratio’ defined as 
(e.g.. [Kalkov et ^l2015h : 


r^{k,z) 


(A^(k,^)|^^l>o.5) 

(A2(k,2)|^J<0.s) 


( 6 ) 


where the averages are over angles. If the signal is isotropic 
then r^(fc, z) will be identically zero, otherwise it can be 
positive or negative. The redshift evolution of the anisotropy 
ratio for scales fc = 0.5,1 Mpc“^ is plotted in Figure[3 The 
source model chosen is SO. We find from the top panels 
that the RSD can cause significant anisotropies (r^ 1.5 

for fc ~ 0.5 Mpc“^) for all the three modelfu. The redshift 


For the models B and C, the level of anisotropies arising from 
the RSD depends on the source model too. The anisotropies are 
caused by the correlation of the 5Tb fluctuations with the dark 
matter density field, which is enhanced when small sources are 
included. For 51 and 52 which do not contain haloes smaller than 
~ 10®Mq, we find that |r^| < 1 from RSD effects (in agreement 
with results of Paper I). 
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Model A 
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Figure 9. Evolution of the ‘anisotropy ratio’ as a function of redshift for the source model S3 at two different scales k = 0.5, lMpc“^. 
The three columns from left to right represent models A, B and C respectively. The top (middle) row represents the case where only the 
RSD (light-cone) effect is included, while the bottom most row represents the case where both the effects are included. 


where the is maximum corresponds to prominent dip in 
the evolution of the power spectrum. 

In contrast, the light-cone effect does not cause any sig¬ 
nificant anisotropy on the signal as can be seen from the mid¬ 
dle panels for the scales considered. The anisotropy, when 
both the effects are included, is thus dominated by the con¬ 
tribution from the RSD. We should mention here that it is 
possible that the anisotropy arising from light-cone effect 
behave differently for larger scales fc < 0.1 Mpc“^. In order 
to have sufficient number of modes at such large scales one 
requires boxes of much larger size which, unfortunately, are 
beyond the scope of this paper. 


4 SUMMARY AND DISCUSSION 

The main focus of the paper is to investigate the impact of 
the light-cone effect on the H l 21-cm signal from the cosmic 
dawn. The 21-cm brightness temperature maps are gener¬ 
ated by post-processing a dark matter A-body simulation 
with an 1-D radiative transfer code. In addition to the usual 
stellar-like sources, we have accounted for X-ray emitting 
mini-quasar like sources. The fluctuations in the spin tem¬ 
perature due to inhomogeneous Lya coupling and heating 
of the IGM not only boost the 21-cm power spectrum but 
introduce several peaks and dips in it when plotted as a 
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function of redshift. The boosted power spectrum, together 
with additional peaks and dips along the redshift axis, makes 
it an ideal case for studying the light-cone effect. Here, we 
calculate the spin temperature fluctuations self consistently 
for various possible source models and reionization histories. 
We then study, for the first time, the light-cone effect on the 
21-cm power spectrum arising from the cosmic dawn when 
the signal is dominated by spin temperature fluctuations 
(unlike other previous studies which concentrated mainly 
on the ionization fluctuations). 

The main findings of our work are summarised below: 

(i) We find that our re sults are consistent with previous 
studies llPatta et al.ll201^ 'l for the model which assumes that 
the H I spin temperature to be fully coupled with the IGM 
temperature and very high compared to the CMB temper¬ 
ature throughout the reionization history (model A). The 
light-cone effect is most significant when the universe is 
^ 15% and 80% ionized as the spherical averaged power 
spectrum is increased/suppressed by a factor of ~ 1.5 and 
~ 0.8 respectively. The effect is found to be minimum while 
the universe is ~ 50% ionized. 

(ii) The light-cone effect is much more dramatic in the 
model B where the inhomogeneous heating is taken into ac¬ 
count. We find that the maximum suppression of the large- 
scale power spectrum due to the light-cone effect occurs 
around the two peaks and enhancement occurs around the 
dip between the two peaks when plotted against 2. We notice 
that the light-cone effect can increase the power spectrum 
by a factor of ~ 1.25 and ~ 3 times around redshift 16 and 
11 respectively at large scales (fc ~ 0.05Mpc“^). The sup¬ 
pression of the power spectrum around the peaks around 
redshifts 8 and 15 can be ~ 20 to 25% at large scales. We 
hnd the enhancement/suppression to be higher for models 
where contributions of small sources have been included (S3 
source model). Unlike the previous studies, we find a signif¬ 
icant light-cone effect at small scales (fc ^ 1 Mpc“^) as well 
at high redshifts. 

In addition, when we calculate the Lya coupling self- 
consistently (model C), we find that the power spectrum 
increases very rapidly at the very beginning of cosmic dawn 
and thus the light-cone effect can enhance the 21-cm power 
spectrum by a factor ~ 1.5 at large scales (fc ~ 0.05 Mpc“^). 
At these scales, the three peaks of the power spectrum are 
suppressed by factors of 0.6, 0.75 and 0.8 around redshifts 
16.5,14 and 8.5 respectively. The light-cone effect can en¬ 
hance the power spectrum in the dips by factors of 2 and 
2.5 around redshift 15.5 and 11 respectively. In general, the 
light-cone effect at small scales is found to be smaller than 
the effect at larger scales, except during the cosmic dawn 
(may increase the power spectrum by factor of few at scales 
like fc ~ 1 Mpc“^). 

(iii) Our results can be understood from the fact that 
the light-cone effect is strong in a situation where the power 
spectrum evolve non-linearly with redshift (or comoving dis¬ 
tance) as any linear evolut ion gets cancelled out to a large 
extent llPatta et al.ll201^ i. Non-linear evolution of the 21- 
cm power spectrum is maximum when it takes a turn. For 
example, for the A model it happens twice, first time when 
iHi ~ 0.8 and second time when xm ~ 0.3 — 0.4. These 
are the places where the light-cone effect is very strong. For 
the B type model, one more peak in the power spectrum 


gets added. Because of this non-linear evolution around the 
heating peak, we see that the light-cone effect becomes very 
strong (power spectrum is suppressed substantially) around 
this peak. Similarly Lya coupling inhomogeneities give rise 
to two more episodes of non-linear evolution of the power 
spectrum (one dip and one peak) in the very beginning of 
cosmic dawn where the light-cone effect is very strong. 

(iv) The large scale {k ~ 0.05 Mpc“^) power spectrum 
with light-cone effect included can differ by ~ —100 to 100 
mK^ from the case where light-cone effect is not included. 
The difference increases at small scales (k > 0.5Mpc“^) to 
the range ~ —250 to 100 mK^. In principle, such a strong ef¬ 
fect should easily be detected by future experiments like the 
SKA. We also notice that the peaks and dips in the power 
spectrum are somewhat smoothe d out due to the light -cone 
effect. It has been suggested fe.g.. lMesinger et al.ll20l3 l that 
the peaks can be used to extract source properties. X-ray 
and Lya background etc. We argue that the light-cone ef¬ 
fect should be considered while extracting those parameters 
from the peak and dip measurements. 

(v) The light-cone effect can, in principle, introduce 
anisotropies in the power spectrum, similar to the RSD 
effects. It is somewhat difficult to predict the large scale 
k < 0.1 Mpc“^ anisotropies because of the lack of num¬ 
ber of modes in the simulation box. However, for interme¬ 
diate scales k ~ 0.5 Mpc“^, we do not find any significant 
anisotropy arising from the light-cone effect. 

The light-cone effect, during the cosmic dawn (ihii 
0.2) is highly dependent on nature of the X-ray sources 
present during that time in the universe. For example, the 
heating of the IGM will be m uch more homogene ous if it is 
mainly driven by the HMXBs llFialkov et al.l[20l3 l. In such a 
scenario the evolution of the signal will be solely dominated 
by the H l neutral fraction and thus the light-cone effect 
during the cosmic dawn will not be as strong as we see here. 
Although, we expect similar light-cone effect if supernova or 
hot inter-stellar gas dominate the X-ray budget since they 
are expected to have similar kind of SED like mini-quasars 
which emit large number of soft X-ray photons. 

Finally, further investigation is required in order to un¬ 
derstand the light-cone effect for different source models and 
reionization scenarios. The main result from this study is 
that the light-cone effect can change the amplitude of the 21- 
cm power spectrum as well as the shape of the evolution of 
the power spectrum at large scale quite substantially. There¬ 
fore, it is important that the effect is incorporated while 
modelling the 21-cm power spectrum, particularly during 
the cosmic dawn and early stages of reionization. 
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